TOMOGRAPHY OF SOLITONS * 

Sergio De Nicola, 1 , Renato Fedele, 2 Margarita A. Man'ko 3 and Vladimir I. Man'ko 3 

1 Istituto di Cibernetica "Eduardo Caianiello" del CNR Comprensorio "A. Olivetti", Fabbr. 70 
Via Campi Flegrei, 34, 1-80078 Pozzuoli (NA ), Italy 
E-mail: s. denicola@cib. na.cnr.it 

■ 2 Dipartimento di Scienze Fisiche, Universita Federico II and INFN, 

o : 

■ Complesso Universitario di M.S. Angelo, 

wo: 

2 ■ Via Cintia, 1-80126 Napoli, Italy 

< 



C3 



Email: renato.fedele@na.infn.it 



\ 3 P. N. Lebedev Physical Institute, 

119991 Moscow, Russia 



Email(s): mmanko@sci.lebedev.ru, manko@sci.lebedev.ru 



CN ' Abstract 
> 

■ We develope the tomographic representation of wavefunctions which are solutions of the gener- 
O , 

alized nonlinear Schrodinger equation (NLSE) and show its connection with the Weyl-Wigner 
map. The generalized NLSE is presented in the form of a nonlinear Fokker-Planck-type equa- 



tion for the standard probability distribution function (marginal distribution). In particular, 
this theory is applied to the envelope solitons, where tomograms for envelope bright solitons of 



^ ■ a wide family of modified NLSE are presented and numerically evalueted. 
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I. INTRODUCTION 



In quantum mechanics and quantum optics as well as in signal analysis, the properties of 
wave functions satisfying the linear Schrodinger equation and the properties of complex analytic 
signals were studied using the tomographic probability distribution functions called tomograms. 
In the tomographic representation, the complex wave function is associated with the standard 
probability distribution by means of invertible map. In signal processing, this method was 
called noncommutative tomography of analytic signal. 

Several tomographic schemes are known, i.e., optical tomography, symplectic tomography, 
photon-number tomography, and spin tomography. An advantage to use the tomographic 
representation for quantum states is related to the possibility of measuring the states, described 
in the standard approach by complex wave functions. The tomographic map of the measured 
tomogram provides the possibility to reconstruct the wave function (up to the constant phase 
factor) of the quantum state. 

Till now the tomographic representation was used for discussing the solutions of linear 
Schrodinger equation of quantum mechanics. 

On the other hand, there exist nonlinear processes described by well-known nonlinear equa- 
tions like nonlinear Schrodinger equation with cubic nonlinearity or more complicated nonlin- 
earities. There exist a huge literature to construct solutions of nonlinear equations and different 
methods to approach the problem. Among solutions of nonlinear equations, there is a set of 
specific soliton solutions. 

Up to our knowledge, till now the tomographic representation was not used to describe the 
soliton solutions of nonlinear equations. We introduce here the tomographic representation of 
the soliton solutions of nonlinear equations and concentrate on tomograms of envelope solitons 
of nonlinear Schrodinger equation. 
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II. WEYL VILLE WIGNER MAP 



Now we review the approach called symplectic tomography of quantum state. The real 
meaning of this scheme is the map of the complex wave function ip(x) of a real variable x 
(—00 < x < 00) onto a family of probability distributions w(X,fi, v) of a random real variable 
X (—00 < X < 00) labeled by two real parameters /i (—00 < fi < 00) and v (—00 < v < 00). 

We use dimensionless variables. The map may be realized by the following steps. First, one 
constructs the density matrix |IJ which is a complex function of two variables x and x' 

p,p(x,x') = ip(x)tfj*(x'). (1) 

Then one uses Wigner- Weyl map of density matrix onto real Wigner function on phase space 
W(q,p) [T2J of two real variables p and q 

W^q,p)= [ Pqp (q + u/2,q-u/2)e- ilm du. (2) 

The Wigner function takes real values. If the wave function is normalized 

J \i)(x)\ 2 dx = 1, 

the Wigner function is also normalized 

J Wj,(q,p) dqdp/2-K = 1. 

Ville @ used this map in the analytic signal theory. 

The inverse of the Fourier transform ([|) defining the Wigner function in terms of density 
matrix reads 

VWOO = -^l w * (^y~ A elv[x ~ x,) d P- ( 3 ) 

One can see that the Wigner function determines the complex wave function up to the constant 
factor 

= (I .p) eipxd p- ( 4 ) 
3 



Modulus of the constant factor \ijj(0)\ is determined by the relationship 



I^(0)| 2 = t^ / W^(0,p) dp. 

271 J 



(5) 



We suppose that ip(0) is not equal to zero. Thus, given Wigner function one can reconstruct 
the complex wave function up to constant phase factor. This means that the Wigner function 
contains the same information which the density matrix does. Also this means that the Wigner 
function contains the same information as the wave function ip(x) does (up to the constant 
phase factor). The Wigner function can be identified with so-called Weyl symbol of density 
operator describing the quantum state. There exist diferent kinds of symbols of operators. As 
it was shown recently [|J tomograms can be also identified with a specific symbol of a density 
operator. 



Let us now construct the tomographic map. To do this, we use the integral Radon transform 
of the Wigner function 



In this formula, the Dirac delta-function term 5(X — fiq — up) collects values of the Wigner 
function W^(q,p) from the line in the phase space which is described by the expression obtained 
by equating the argument of Dirac delta-function to zero. One can prove || that for normal- 
ized Wigner function the function w(X,/a, v) is normalized probability distribution function of 
random variable X (called tomogram of the Wigner function), i.e., one has 



III. TOMOGRAPHIC MAP 




(6) 




(7) 



This tomogram was called symplectic tomogram ||, since it is related to linear symplectic 



transform in the phase space. 



There is another form of the map (^) given in terms of Fourier transform (we used Fourier 
transform of delta-term in (|6|)) 
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/elk dc£ dp 
W^{q,p) exp [ik(X -nq- up)] . (8) 

The Fourier integral form of the tomogram (§) gives the possibility to get easily the inverse 
transform 

f dX dfidu 
W{q,p) = J w{X,fi, v) exp [i(X - fiq - up)] — . (9) 

Thus, given the symplectic tomogram w(X,fi,u) one can find the Wigner function W(q,p). 
Using known relationships ([[]), (|2]) and (Q) of the Wigner function W(q,p) and the wave function 
i/j(x), one can obtain the expression for tomograms in terms of the wave function M 

'%H 2 .Xy" 



w(X, /I, u) 



2ttM 



j tp(y) exp y 2 - i dy 



(10) 



This formula shows that for a complex function i/j(x) one can find the tomogram. Formula (|10|) 
was found for the quantum wave function and for noncommutative tomography of the analytic 
signal || . But it can be used for other arbitrary aims as well. The goal which we are going to 
reach is to use this formula for the description of solitons. 

IV. TOMOGRAM AS FOURIER TRANSFORM OF A CHIRPED SOLITON 

There are some properties of tomograms to be used. The homogeneity property follows 
from relations (||) and ( JTOj ) 

w(XX, Xfi, = -p-w(X, fi, v). (11) 

I'M 

This means that, in reality, the tomogram is the function of two real variables. For example, 
one can take 

/i = cos 9, u = sm6. (12) 

In this case, the symplectic tomogram is given by the formula 

w(X,9)= [ W(q,p)5(X - qcos6 ~psm6)^^-. (13) 
J 2n 



This relation of the Wigner function W(q,p) and the tomographic probability w(X,9) (called 
tomogram of the optical tomography scheme) was used in the signal theory in ||10|| . The relation 
of the optical tomogram w(X, 9) to the wave function reads 



w(X,9) 



2tt \ sin 9\ 



J i/>(y) exp Q cot 9y 2 - dy 



(14) 



The integrand in the above expression (|lj) is similar to the Green function of the quantum 
harmonic oscillator and the tomogram coincides with modulus squared of fractional Fourier 
transform of ip{y) WM - The tomogram w(X, 9) is used in quantum optics in the scheme of mea- 
suring quantum photon states by means of the so-called optical homodyne tomography ||12|| . It 
was also discussed in the context of quantum-optics measurements in [|13|]. Thus the tomogram 
w(X,fi,u) determines completely the Wigner function and the complex function ip(x) (up to 
the phase factor), if it is known for real parameters satisfying the constraint 



fi 2 + v 2 = 1. 



(15) 



But one can use other constraints too. Thus, the homogeneity property ( |TT| ) implies that 
the particular values of the symplectic tomogram, e.g., w(l,fi, is), w(X, 1, v) and w(X, /i, 1) 
determine the whole tomogram and, consequently, the complex function ip(x) (up to the phase 
factor) and the Wigner function W(q,p) completely. 

Another property can be obtained by change of variables in (|10|) 



v 



(16) 



which gives the following expression for the tomogram 



w(X, fi, v) 



2tt 



(17) 



where the function if)(z,fi,v) is the function describing "the chirped soliton" (we mean that 
ip(y) in ( |i~0| ) is considered as a soliton solution of a nonlinear equation) 



il>(z, n, v) = %l){zv) exp ( - fiuz 2 



Expression (|17|) is convenient for numerical calculations because it gives the tomogram in 
terms of the standard Fourier transform of the chirped soliton. In terms of optical tomogram, 
the expression can be rewritten as 



w(X,6) 



sine 
~27 



iP(z,6)e- iXz dz 



(19) 



where the chirped soliton has the form 



iz 



tp(z, 6) = ip(z sin 6) exp I — sin 26 ) . 



(20) 



The formula obtained is used below to make plots of solitons. 



V. NONLINEAR EQUATIONS IN TOMOGRAPHIC AND 
WEYL WIGNER MOYAL REPRESENTATIONS 



Let us consider the following generalized nonlinear Schrodinger equation (NLSE): 



. dip 1 d 2 ip 



ds 



2 dx 2 



+ U 



*P, 



(21) 



where s and x are the time-like and space-like variables and ip = ip(x, s) is a complex wave func- 
tion describing the system's evolution in the configuration space; U = U [l^l 2 ] is an arbitrary 
real functional of \ip\ 2 - 

In this section, we derive the evolution equations for solitons in the phase-space represen- 
tation. For the case of the cubic NLSE 



.dip 1 d 2 ip 



ds 2 dx 2 

the density matrix (JJ) satisfies the following evolution equation: 



(22) 



.dp(x, x', s) 



1 / d 2 



d 2 



ds 



2 \ dx 2 dx' 2 



p(x,x', s) + 



5(x - y)p(x,y,s) dy 



V 



5{x' - y)p(x',y,s) dy 



p(x,x',s) 



(23) 



where 



V 



S(x-y)p(x,y,s) dy 



J S(x-y)p(x,y,s)dy 



(24) 



is the potential-energy functional of the density matrix for nonlinear Schrodinger equation 
with cubic nonlinearity. 

The transition to the evolution equation for the Wigner function can be done using the 
standard algebra, which provides the following recepie. One has to make in (|23|) the replacement 
p — > W along with the following replacement: 







p(x, x') 



dx 

x p (x, x') 
x' p (x, x') 



ld_ 

2dq 



ip) W (q, p) 



(25) 



It provides the following Moyal-like form of the nonlinear equation (|23|) for the Wigner function 
dW(q,p, s) 



ds 



dW(q,p } s) 1 j 
-p ^ + - < V 

oq ! 



/ id id 
T + 2V g+ 2^' S , 



c.c. \W(q,p,s). (26) 



Here the arguments of the potential energy are replaced by the operators which act onto the 
Wigner function. For the cubic NLSE 



V(z) = q z, 



(27) 



one has a simple form of the equation 

dW(q,p,s) d TTr , . _ / id id \ rTr . . 
ds + P dq ^ P ' S) ~ 2qolmp [ q+ 2dp~ ,q+ 2dp~ ,S ) ^' P ' S) = 



Using the relation 



p{x,x) = J W(x,p) y 



one has 



dW(q, p, s) dW(q, p, s) 
■ +p- 



r I id \ dP 
2g Im j w(q + -—,P,s\—W(q,p,s) = 0. 



ds 



dq 



' dp 



(28) 



(29) 



(30) 



In (|26|), (0), and ( j30| ) the arguments of the density matrix and Wigner function are replaced 
by operators and the operators act on the Wigner function itself. Equations fl26|) and (|28|) can 



be presented in the form of Moyal-like series ]14[]. Equation fl26|) for an arbitrary nonlinear 
potential V(z) can be written in terms of the functional partial differential equation for the 
Wigner function only 

Via r, S \ ( T r ( ■if) \ JP~\ 1 

W(q,p,s) = 0. (31) 



dW(q, P,s)_ +p dW(q,p, s) _ J v 



rxrf i d „ \ dP 



ds ' dq 

Now we consider the relation of the density matrix, Wigner function and tomogram 
W(q, p) = — [ w(X, p, v) exp [—i(pq + vp — X)] dp, dv dX ; 

Z7T J 



p(X, X') = ±-l w(Y, p,X- X') exp 

Z7T J 



i[Y-p 



X + X r 



dp dY. 



(32) 
(33) 



In view of this relation, the following rule for substitution in the evolution equation 



p(x, x' , s) — > w(X, p, v, s) 
\dXi 



xp 



i d 
dp + 2 U dX 



d_ 

dx^ 



x p 




p d 
2dX^ 1 [dX 



d \ 1 d 



dX 



-1 



d 



dv 




w: 



UK 



(34) 



i 

'dp ~2 U dX 



dx 



p d 

2dX +t [dX 



d \ 1 d_ 
dv 



w; 



w 



provides the tomographic form of the nonlinear equation under consideration 



dw(X, p,v, s) dw(X, p,v, s) 



ds ' dv 

— 2 Im V < J w(y, p', 0, s) exp 



i\y + P 



'j^y 1 d_ 

.dXj 



i d 
dp ~2 V ~dX 



dy dp' 
2tt 



w(X, p, v, s) — 0. (35) 



In the above equation (|35|) , one has the integro-differential operator in the exponent, which 
acts on the tomographic probability function. The integral operator (d/dX)' 1 is defined by 
the following action on Fourier component of the function f{X): 



e ikx dk. 

ik 



For the case of the cubic NLSE, one has 



dw(X, fj,,u,s) dw(X,n,v,s) 



ds ' dv 

-2q Im < J w(y, //, 0, s) exp 



ily + fi 



'_d_ 

dX 



-1 



did 
djl ~2 U dX 



dy djx' 
2it 



w(X, fi, v, s) = 0. (36) 



It should be pointed out that Eq. fl5E|) has the solutions which in the course of the evolution 
process preserve the positivity and normalization. Thus the soliton solutions of the nonlinear 
equations can be mapped onto probability distribution functions. 

The meaning of the probability distributions is the following. If x is a coordinate and p 
is the momentum, the value X = fix + up is the position in the reference frame in the phase 
space (x,p), the reference frame being scaled and rotated. The parameters of the scaling A and 
rotation 6 are determined by the real parameters \x and v, namely, fi = e x cos 9 and v = e~ x sin 9. 

The probability distributions w(X, fi, v) determine soliton solutions ip{x) in the correspond- 
ing representation. Consequently, in the tomographic representation soliton solutions of nonlin- 
ear dynamic systems are solutions of generalized Fokker-Planck-type equations for the standard 
probability distributions. Such representation can be useful from mathematical point of view 
since analysis of probabilities and their assymptotics can be additionally incorporated (using 
existing theorems) on the behaviour of the probability distribution functions. 

VI. EXAMPLES 



In this section, we study bright soliton in both the tomographic and Weyl-Wigner repre- 
sentations. The envelope bright soliton of the cubic NLSE, i.e. 

(2\E\ 



00 



v 1/2 








sech 




exp 


i 



V x -(E + V*/2) s 



(37) 



where E is a negative real constant, Vq is an arbitrary real constant and £ = x — VqS (see f.i., 
Ref. ||15| ). Thus, the corresponding optical tomogram is given by the formula 

IS sin #1 



w h (X,9,s) 



0o F 



sech 



2\E\(ysm0-V o s) 



x 



10 



x exp < i 



V y sin 9 



% sin 28 
+ — y 



2 iXy dy 



(38) 



It has been recently shown that the following modified NLSE ([/[|-?/>| 2 ] = go I 



ds 



d 2 ^ 
2~dx 2 ~ 



(39) 



for go < and any real positive value of j3, has the following envelope soliton-like solutions 

H 1/2/3 

exp 



s) 



sech 1//3 



(40) 



where the real numbers Vo and E are arbitrary and negative, respectively, and still £ = x — Vos 
(note that Vo is the soliton velocity). It should be noted that the case (3 = 1 (ordinary envelope 
bright soliton of the cubic NLSE |17| ) can be very easily recovered [|T6| . 
The Wigner function of bright soliton for Vo = is given by the formula 



\E\ f 
W(x,p) = - — - / sech 

I go J 



2\E\(x + ^ 



sech 



a 



- / i x - - 



exp~ ip " du. 



(41) 



The tomogram of the soliton solution of generalized nonlinear Schrodinger equation is given by 



w(X, /i, v) 



= 0) 




1 


l\E\(l + /3) 


2tx\v\ 


\Qo\ 



sec 



N2\e\ y 



ifi 2 iXy\ 

ex P 1 7rv d y 



2u 



V ) 



(42) 



The Wigner function of the soliton solution of generalized nonlinear Schrodinger equation is 
given by the formula (Vo = 0) 



W(x,p) 



\E\(1 + P) 
\lo\ 



V/3 



sech 1//3 



P\/2\E\(x + - 



u 



sec. 



/3y/2\E\[x-- 



a 



exp-' ipu du. 
(43) 



VII. SOME NUMERICAL COMPUTATIONS 



3D Plots and density plots of both tomograms and Wigner functions for bright solitons with 
(3 = 0.5, 1.0, 2.0 and 2.5 are given in Fig.l - Fig.5. We have fixed the free parameters as follows: 
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Vo = 0, E = — 1, and go = — 1- Fig.l represents the 3D plots of the tomogram of the solitons 
for the different values of (3. The corresponding density plots are displayed in Fig. 2 . Fig. 3 
displays the 3D plots of the Wigner quasidistribution of the solitons for the diverse values of j3. 
For such solitons, the phase-space regions where the Wigner function is negative are not easily 
visible in Fig.3, whilst the corresponding density plots displayed by Fig.4 show clearly this 
behaviour. They correspond to the black and grey regions. The "deepness" of the negativity 
is represented in scale of grey. Black regions correspond to the deepest negative parts. The 
size of these negative parts of the Wigner function are, for instance, clearly represented by the 
cross section shown in Fig. 5 at p = 2 for the diverse values of j3. 



We discussed the tomography of solitons considering the tomograms as additional char- 
acteristics of the soliton solutions of nonlinear dynamic equations. But there exists another 
experimental aspect of tomograms because the tomograms can be directly measured in differ- 
ent situations. Thus the problem of tomography of some phenomena which is described by a 
complex function ip(x) is equivalent to the problem of measuring the amplitude 1^(^)1 an d the 
phase <^(x) of the complex function ip(x) = \ip(x)\ expi<p(x). The function ip(x) can describe a 
soliton but also this function can describe some signal connected with different processes, e.g., 
in optical fibers, in plasma, etc. In all processes where one needs to measure the amplitude 
and phase by measuring experimentally only intensities, the tomography can be used as an 
instrument for achieving this aim. We consider two different possibilities which are based on 
the symplectic tomogram flT0|). The density matrix can be reconstructed either by measuring 
the tomogram of the optical tomography scheme w(X, 6) or tomogram w(X, l,v). 

Now we show that both tomograms can be obtained in two different and realizable processes. 
The optical tomogram can be rewritten in terms of fractional Fourier transform (which is 
reduced to the Green function of the harmonic oscillator) In fact, one has 



VIII. REMARKS AND PERSPECTIVES 



1 




iXy] 



2 



w(X,G) 



\J 2ni sin 



sin 



dy 



(44) 
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In this formula, we take v = sinO and /i = cosG. The phase factor expiX 2 /2cotG does not 
change the value of the tomogram. 

The tomogram presented in such form coincides with the value of the wave function at the 
point x at the time moment t if the initial value of the wave function at the time moment t = 
is equal to ipiu)- This means that to reconstruct the initial value of the wave function ijj(x) 
including both the amplitude ^(x)! and the phase <p(x) 



one can measure the tomogram, i.e., amplitude squared of the wave function which evolves 
in the quadratic potential well. This situation can be perfectly done for optical fibers with 
a parabolic profile of the refractive index called "selfoc" (linear propagation). In fact, the 
light beams in optical fibers obey to the Schrodinger-like equation which follows from the 
Helmholtz equation in the Fock-Leontovich approximation |Tj|. But time t in the Schrodinger 
equation is replaced by the longitudinal coordinate z and the Planck's constant is replaced 
by the wavelength. Thus, to measure the input field amplitude and phase, it is sufficient to 
measure the tomogram which is the field intensity in each cross-section of the fiber given by 
longitudinal coordinate < z < 2tt. 

Another possibility is related to the formula 



This formula is equivalent to formula (|T0|) in which we took [i = 1 and added unessential phase 



due to homogeneity the tomogram w(X, 1, v) is equivalent to the tomogram w(X,fi, u), while 
measuring the intensity of free propagating signal one measures both the phase and amplitude of 
the input signal ip{y). If one measures the field in optical fiber, the structure of the output field 
can be evaluated by measuring the free propagation of the beam. There is a pecularity in using 
formula (f45l). For complete reconstructing the amplitude, one needs to know the intensity for 



ip(x) = \ip(x)\ expiip(x), 




2 




(45) 



factor exp Thus the tomogram for "time moment" v is equivalent to the intensity of the 
free propagating signal. In fact, the kernel in (f45D is the Green function of a free particle. Since 
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arbitrary large values of time (or longitudinal coordinate z). Practically the length or duration 
can be chosen to fit appropriate accuracy of the measurement. As an example, Fig. 6 shows the 
tomographic map and the corresponding density plot for the case of free propagation of a pulse 
whose initial profile is solitonlike with (5 — 1.0 (also here we have fixed the free parameters as 
follows: V = 0, E = —1, and q = —1). 

IX. CONCLUSIONS 

We introduced the tomographic probability distribution associated with soliton solutions of 
nonlinear equations. 

Nonlinear dynamical equations like nonlinear Schrodinger equation were presented in the 
form of equation (a nonlinear generalization of the Fokker-Planck equation) for the standard 
probability distribution function. 

Specific cases of solitons for a wide family of modified NLSE were studied in the tomographic 
representation explicitly. 

The possibility to use tomograms to reconstruct the phase of linear or nonlinear signals by 
measuring the signals' intensities was finally discussed. 
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